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X-ray observations of cluster outskirts: 
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Abstract. 

Past and current X-ray mission allow us to observe only a fraction of the volume occupied 
by the ICM. After reviewing the state of the art of cluster outskirts observations we discuss 
some important constraints that should be met when designing an experiment to measure X- 
ray emission out to the virial radius. From what we can surmise WFXT is already designed 
to meet most of the requirements and should have no major difficulty in accommodating the 
remaining few. 

Key words, galaxies: cluster: general - galaxies: fundamental parameters - intergalactic 
medium - X-ray: galaxies - cosmology: observations - dark matter 
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1. Introduction 

Galaxy clusters form through the hierarchi- 
cal accretion of cosmic matter The end prod- 
ucts of this process are virialized structures 
that feature, in the X-ray band, similar ra- 
dial profiles of the surface brightness Sb (e.g. 
Vikhlinin et al. 1999, Neumann 2005, Ettori 
& Balestra 2009) and of the plasma temper- 
ature Tgas (e.g. Allen et al. 2001, Vikhlinin 
et al. 2005, Leccardi & Molendi 2008). Such 
measurements have definitely improved in re- 
cent years thanks to the arcsec resolution and 
large collecting area of the present X-ray satel- 
lites, like Chandra and XMM-Newton, but still 
remain difficult, in particular in the outskirts, 
because of the low surface brightness asso- 
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ciated to these regions. Present observations 
provide routinely reasonable estimates of the 
gas density, Wg^s, and temperature, Tgas, up to 
about 7^2500 (~ 0.3/? 200; R^ is defined as the ra- 
dius of the sphere that encloses a mean mass 
density of A times the critical density at the 
cluster's redshift; /?200 defines approximately 
the virialized region in galaxy clusters). Only 
few cases provide meaningful measurements at 
R5QQ {~ 0.7/?20o) and beyond (e.g. Vikhlinin et 
al. 2005, Leccardi & Molendi 2008, Neumann 
2005, Ettori & Balestra 2009). Consequently, 
more than two-thirds of the typical cluster vol- 
ume, just where primordial gas is accreting and 
dark matter halo is forming, is still unknown 
for what concerns both its mass distribution 
and its thermodynamical properties. This poses 
a significant limitation in our ability to charac- 
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terize the physical processes presiding over the 
formation and evolution of clusters and to use 
clusters as cosmological tools, as also outlined 
in the Scientific Justification for the WFXT 
(Giacconi et al. 2009). Indeed the character- 
ization of thermodynamic properties at large 
radii would allow us to provide constraints on 
the virialization process, while measures of the 
metal abundance would allow us to gain insight 
on the enrichment processes occurring in clus- 
ters (e.g. Fabjan et al. 2010). Morever the X- 
ray emission at large radii could also be used 
to improve significantly measures of the gas 
and total gravitating masses thereby opening 
the way to a more accurate use of galaxy clus- 
ters as cosmological probes (e.g. Volt 2005). 

In these proceedings, we take stock of the 
situation on cluster outskirts and suggest how 
to make progress. In Sect. |2] we provide an 
observational overview of cuiTently available 
measures of cluster outer regions, while in 
Sect. [3] we discuss some important constraints 
that should be met when designing an exper- 
iment to measure X-ray emission out to the 
virial radius. In Sect.|4] we present an overview 
of future missions which have cluster outskirts 
observations as one of their goals, our main re- 
sults are recapitulated in Sect. |5] 

A Hubble constant of 70 hjo km s"' Mpc"' 
in a flat universe with Q,n equals to 0.3 is as- 
sumed throughout this manuscript. 



2. What we know of cluster outskirts 

2.1. Surface brightness and gas density 
profiles 

The X-ray surface brightness is a quantity 
much easier to characterize than the temper- 
ature and it is still rich in physical informa- 
tion being proportional to the emission mea- 
sure, i.e. to the gas density, of the emitting 
source. Recent work focused on a few local 
bright objects for which ROSAT PSPC obser- 
vations with low cosmic background and large 
field of view have allowed to recover the X- 
ray surface brightness profile over a significant 
fraction of the virial radius (Vikhlinin et al. 
1999, Neumann 2005). 



In Ettori & Balestra (2009), we study the 
surface brightness profiles extracted from a 
sample of hot (Tgas > 3 keV), high-redshift 
(0.3 < z < 1.3) galaxy clusters observed 
with Chandra and described in Balestra et al. 
(2007). A local background, B, was defined 
for each exposure by considering a region far 
from the X-ray center that covered a signifi- 
cant portion of the exposed CCD with negli- 
gible cluster emission. We define the "signal- 
to-noise" ratio, S2N, to be the ratio of the ob- 
served surface brightness value in each radial 
bin, S bir), after subtraction of the estimated 
background, B, to the Poissonian error in the 
evaluated surface brightness, ei(r), summed in 
quadrature with the error in the background, 

6b: SlNir) = [Stir) - B] / ^etir)^ + ej. The 

outer radius at which the signal-to-noise ratio 
remained above 2 was defined to be the limit of 
the extension of the detectable X-ray emission, 
Rs2N- We estimated i? 200 using both a^S-model 
that reproduces the surface brightness profiles 
and the scaling relation quoted in eq. [T] and 
selected the 11 objects with Rsin/Rkk) > 0.7 
to investigate the X-ray surface-brightness pro- 
files of massive clusters at r > R;,oo ~ 0.7/?2oo- 
Examples of the analyzed dataset are shown 
in Fig. [1] We performed a linear least-squares 
fit between the logarithmic values of the ra- 
dial bins and the background-subtracted X-ray 
surface brightness. Overall, the error-weighted 
mean slope is -2.91 (with a standard deviation 
in the distribution of 0.46) at r > 0.2/? 200 and 
-3.59(0.75) at r > 0.4/?20o- For the only 3 ob- 
jects for which a fit between 0.5/?200 and Rsin 
was possible, we measured a further steep- 
ening of the profiles, with a mean slope of 
-4.43 and a standard deviation of 0.83. We also 
fitted linearly the derivative of the logarithm 
S bir) over the radial range 0.1/? 200 - Rsin, ex- 
cluding in this way the influence of the core 
emission. The average (and standard deviation 
cr) values of the extrapolated slopes are then 
-3.15(0.46), -3.86(0.70), and -4.31(0.87) at 
0.4/?2oo, 0.7/?2oo and /?200, respectively. 

These values are comparable to what has 
been obtained in recent analyses. Vikhlinin et 
al. (1999) find that a /3-model with/3 = 0.65 - 
0.85 describes the surface brightness profiles 
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Fig. 1. From left to right: Example of a surface brightness profile with the fitted background (hor- 
izontal dotted line) and the radius R200 (vertical dashed line); the signal-to-noise profile evaluated as 
52A' = (Sj, - B)le, where the error e is the sum in quadrature of the Poissonian error in the radial counts 
and the uncertainties in the fitted background, B\ the best-fit values of the slope of the surface brightness 
profile as a function of r/R2oo- These values are estimated over 6 radial bins (thick horizontal solid line: 
the slope evaluated between 0.4 x R20Q and Rsin with a minimum of 3 radial bins; dashed line: best-fit of 
dln(Si,)/dln(r/R2m) with the functional form sq + si ln(r//?20o) over the radial range 0.1 xR2m-Rs2N, with 
the best-fit parameters quoted in Table 3 of Ettori & Balestra 2009). 



in the range 0.3 - IRm) of 39 massive local 
galaxy clusters observed with ROSAT PSPC. 
For ay6-model with x = r/r^, dlnS b/dlnx - 
(1 - 6f^)x^/(l + x^) and 51n«gas/51n;ic - 
-3/3x2/(1 + x^), impling that/3 = 0.65 - 0.85 
corresponds to a logarithmic slope of the sur- 
face brightness of -2.9/ - 4.1, that is a range 
that includes our estimates. Neumann (2005) 
finds that the stacked profiles of few massive 
nearby systems located in regions at low (< 
6 X 10^" cm"2) Galactic absorption observed 
with ROSAT PSPC still provide values of /3 
around 0.8 at R200, with a power-law slope that 
increases from -3 when the fit is performed 
over the radial range [0.1, l]/?200 to -5.7;^|2 
over [0.7,1. 2]/?2()0. 

These observational results are supported 
from the hydrodynamical simulations of X- 
ray emitting galaxy clusters performed with 
the Tree-i-SPH code GADGET-2 (Roncarelli et 
al. 2006; see, e.g.. Fig. |2]i. In the most mas- 
sive systems, we measured a steepening of 
Shir), independently from the physics adopted 
to treat the baryonic component, with a slope 
of -4, -4.5, -5.2 when estimated in the radial 
range 0.3 - 1 .2/?2()o, 0.7-1 .2/?20(), 1 -2 - 2.7;?200, 
respectively. In particular, we note the good 
agreement between the slope of the simulated 



surface brightness profile of the representative 
massive cluster in the radial bin 0.7 - 1.2/? 200 
(see values of bA in Table 4 of Roncarelli et al. 
2006 ranging between -4.29 and -4.54) and 
the mean extrapolated value at /?2(K) of -4.43 
measured in the Chandra dataset. 



2.2. Temperature and metallicity profiles 

Early attempts to produce temperature pro- 
files were made with the ROSAT PSPC, these 
were mostly limited to low mass systems 
(e.g. David et al. 1996) where the tempera- 
tures were within reach of the PSPC soft re- 
sponse. Resolved spectroscopy of hot systems 
began with the coming into operation of ASCA 
(1994) and BeppoSAX (1996). Both missions 
enjoyed a relatively low instrumental back- 
ground, which was a considerable asset when 
extending measures out to large radii, however 
they both suffered from limited spatially res- 
olution. The situation was somewhat less se- 
vere with the BeppoSAX MFCS than with the 
ASCA GIS since the former had a factor of 
2 better angular resolution and a modest en- 
ergy dependence in the PSF These difficulties 
led to substantial differences in temperature 
measures, while on the one side Markevitch 
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Fig. 2. From left to right; upper panels: Maps of the projected gas density, mass-weighted temperature 
and soft (0.5-2 keV) X-ray emission. The circles indicate the virial radius. The size of the side of each 
map is 12 Mpc, so they cover roughly up to 2.5i?200- Bottom panels: Comparison between the gas density, 
mass-weighted temperature, soft X-ray surface brightness profiles for a cluster with My„ 2 x IO'^Mq sim- 
ulated by using 4 different physical models (Roncarelli et al. 2006). A dashed line indicates the functional 
from Vikhlinin et al. (2006, eq. 9) that well reproduces the behavior of the temperature profile of nearby 
bright galaxy clusters observed with Chandra. The extragalactic unresolved background from Hickox & 
Markevitch (2006) in the soft X-ray band is indicated by the shaded region. 



et al. (1998) using ASCA and De Grandi & 
Molendi (2002) using BeppoSAX MECS found 
evidence of declining temperature profiles, on 
the other. White (2000) using ASCA and Irwin 
et al. (1999) using BeppoSAX data found flat 
temperature profiles. The situation was some- 
what clearer on abundance profiles were work- 
ers using ASCA (e.g. Finoguenov et al. 2000) 
and BeppoSAX data (De Grandi & Molendi 
2001) consistently found evidence that cool 
core systems featured more centrally peaked 
profiles than NCC system. The coming into op- 
eration of the second generation of medium en- 
ergy X-ray telescopes, namely XMM-Newton 
and Chandra, both characterized by substan- 
tially better spatial resolution, allowed more 
direct measures of the temperature profiles. 
The new Chandra (Vikhlinin et al. 2005) and 



XMM-Newton measurements (e.g. Pratt et al. 
2007, Snowden et al. 2008) confirmed the pres- 
ence of the temperature gradients measured 
with ASCA and BeppoSAX. In a detailed study 
of a sample of 44 objects observed with XMM- 
Newton (Leccardi & Molendi 2008) we found 
that temperature measurements could be ex- 
tended out to about 0.7/?i8o (see Fig.[3]l. Since 
the major obstacle to the extension of measure- 
ments to large radii was the high background, 
most importantly the instrumental component, 
we adopted the source over background cri- 
terion originally introduced in De Grandi & 
Molendi (2002) to decide where to stop mea- 
suring profiles. The source to background ra- 
tio, defined as 7^, where I^ou and Ii,kg are the 
source and background intensities respectively. 
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Fig. 3. Left panel: mean temperature profiles obtained from Leccardi & Molendi (2008; LM08, black cir- 
cles), De Grandi & Molendi (2002; DM02, blue squares), Vikhlinin et al. (2005; V05, red upward triangles) 
and Pratt et al. (2007; P07, green diamonds). All profiles are rescaled by kTu and Risq. The dashed line 
shows the best fit with a linear model beyond 0.2 Rigo and is drawn to guide the eye. lower panel: residuals 
with respect to the linear model. The LM08 profile is the flattest one. Right panel: mean metallicity profiles 
obtained from Leccardi & Molendi (2008b; LM08, circles), De Grandi et al. (2004; DM04, triangle) and 
Baldi et al. (2007; BA07, squares). Abundances are expressed in Anders & Grevesse (1989) solar values 
and radii in units of RigQ. The radii have been slightly offset in the plot for clarity. 



should not be confused with the signal to noise 



ratio defined as 



t, where t is the ex- 



(/,„„+/hs)"^ 

posure time. While the latter ratio is associated 
to the statistical error and therefore increases 
with exposure time, the former is associated to 
the systematic error and does not depend on the 
exposure time. Through a series of tests (see 
Sect. 5.2.1. and Fig. 11 of Leccardi & Molendi 
2008) we determined that measurements could 
be trusted out to radii where the source to 
background ratio in the 0.7-10 keV band re- 
mained above a threshold of 0.6. In Leccardi & 
Molendi (2008) we made use for the first time 
of extensive simulations to estimate the impact 
of systematic errors on the measurements, part 
of the expertise we have acquired from that 
work has been used to perform the simulations 
discussed in Sect. 13.41 In the left panel of Fig. [3] 
we show a compilation of mean temperature 
profiles from different missions, all show evi- 
dence of a decline of the temperature beyond 
0.2/?igo. Interestingly, as a result of the correc- 
tion for systematic that we applied to our pro- 
file (see Sect. 5.3 and Fig. 14 of Leccardi & 
Molendi 2008) ours is the flattest amongst the 



profiles shown in the left panel of Fig. [3] The 
measurement of the metal abundance profile 
extends to radii that are somewhat smaller than 
those reached by the temperature profiles, this 
is because the most prominent emission line, 
the Fe Ka, is located in the high energy part 
of the spectrum where the instrumental back- 
ground is particularly strong. In the right panel 
of Fig. [3] we show the mean abundance profile 
measured with different satellites. The flatten- 
ing of the profiles beyond 0.27? i go is most likely 
indicative of an early enrichment of the ICM 
(Fabjan et al. 2010). 

Unfortunately the high orbit of the XMM- 
Newton and Chandra satellites, as well as 
the fact that the design of the satellites was 
driven by scientific objectives other than the 
characterization of low surface brightness re- 
gions, led to a substantially higher and more 
variable background than with the previous 
satellite generation, thereby limiting the ex- 
ploration of the temperature and abundance 
profiles to roughly the same regions already 
investigated with ASCA and BeppoSAX (see 
Fig.|3]l. Recently measures of temperature pro- 
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files have been made with the Suzaku X- 
ray imaging spectrometer (XIS). Ahhough not 
ideal for cluster measurements, the XIS fea- 
tures a poor PSF and a small FOV, it does enjoy 
the considerable advantage of the modest back- 
ground associated to the low earth orbit. The 
measures have been conducted on a handful of 
systems (A2204, Reiprich et al. 2009; A1795, 
Bautz et al. 2009; PKS0745-191, George et al. 
2009; A1413, Hoshino et al. 2010) and extend 
beyond the regions explored with Chandra and 
XMM-Newton. However, the characterization 
is a limited one at best: only parts of the outer- 
most annuli are explored and both radial bins 
and error bars are large. Moreover there are 
concerns as to the reliability of the measure- 
ments themselves. All measured temperature 
profiles are steeper than those predicted by 
simulations. This is particularly true of A1795 
and PKS0745-191, where the temperature and 
the surface brightness are respectively steeper 
and flatter than those predicted by simulations. 
Consequently entropy profiles are flatter and, 
in the case of PKS0745-191, it features an in- 
version around 0.6/?200, that could be associ- 
ated to the presence of non virialized gas or, 
alternatively, to problems in the characteriza- 
tion of the source spectrum. 

3. How we can map out to /?2oo 



From the discussion in Sect. 12.21 it is rather 
obvious that past X-ray mission were not op- 
timized for the spectral characterization of 
the low surface brightness emission typical of 
cluster outer regions. In this section we dis- 
cuss how to design an experiment character- 
ized by high sensitivity to low surface bright- 
ness emission. The sensitivity depends upon: 
1) the surface brightness of the source, S h, that 
scales with effective area of the experiment, 
Ae\ 2) the solid angle covered by the field of 
view (FOV), Q; 3) the surface brightness of the 
background, B. The quantity that needs to be 
maximized is then: 

27r]^®"'°'A£(6')0£/6l 

s , 

B 

where 6 is the ofi'-axis angle and the inte- 
gration is extended over the full FOV, i.e. 



2ti J 0d9 - Q.. Therefore one needs to maxi- 

mize the numerator, 2n J AE(0)6d6, a quan- 
tity that is often refeiTed to as "grasp", and 
minimize the background^. To go well beyond 
what has been achieved with the instrumenta- 
tion that has been designed thus far one needs 
to operate at three different levels: 1) the exper- 
iment design; 2) the observational strategy; 3) 
the data analysis strategy. 

3.1. Experiment design 

Let us start by considering the background and 
in particular the instrumental background, i.e. 
the part of the background that is not associated 
to genuine cosmic X-ray photons. A few things 
can be easily infeiTed by comparing back- 
ground spectra from different mission. In Fig.|4] 
we report a recent compilation of such spectra 
from Hah et al. (2008). We note that: 1) front 
illuminated CCDs have lower background than 
background illuminated ones and that 2) the 
background on the low earth orbit is smaller 
than that in the high orbit. In this respect it 
is particularly instructive to compare the EPIC 
MOS with the SWIFT XRT background, since 
we are dealing with virtually the same detec- 
tor in a high and low earth orbit. As shown 
in Hall et al. (2008), the SWIFT XRT back- 
ground is about a factor 3 lower than the EPIC 
MOS background. Thus, from the inspection 
of Fig. |4] we learn that to keep the instrumental 
background low it is preferable to employ front 
illuminated CCDs on a low earth orbit. There 
are other issues that should be kept in mind: 1) 
a non-negligible fraction (say 15%) of the de- 
tector should be shielded from the sky, this will 
allow to constantly monitor the intensity of the 
instrumental background; 2) a tilted CCD con- 
figuration which allows to improve the imag- 
ing, will result in fluorescence Si line emission 
inhomogeneous distributed on the FOV, some- 
thing similar is observed on MOS EPIC, this 
can be minimized by studying the most appro- 
priate configuration; 3) while active shielding 



' A substantial fraction of the background is of 
instrumental origin. This part and Ae scale with the 
square of the focal length. Both of them appears in 
the quantity s that has to be maximized 
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Fig. 4. Left panel: compilation of instrumental background spectra, only continuum components, for var- 
ious X-ray missions equipped with CCD detectors from Hall et al. (2008). Right panel: Resolved fraction 
of extragalactic cosmic X-ray background in the 0.5-2 keV band as a function of angular resolution.The 
total background intensity is derived from measures by De Luca & Molendi (2004) and Mc Gammon et al. 
(2002). The LogN-logS is taken from Moretti et al. (2003). We also include an euclidian component to ac- 
count for the unresolved 20% of the cosmic X-ray background (CXB). The normalization of this component 
is conservatively chosen in such a way that about 10% of the CXB is found at fluxes that are sufficiently 
small to efficiently mimic a diffuse component. The angular resolution necessary to reach a given flux limit 
is obtained by imposing that the source density at that flux limit is such that there is 1 source every 20 
angular resolution elements, where the angular resolution element is a circle with a radius equal to the 
half-power-radius of the PSF. 



cannot be applied as long as the detector is 
a CCD, passive shielding can and should be 
considered. Most importantly the whole instru- 
mental background issue should be addressed 
from a global point of view. Detailed simula- 
tions of the physical interaction between par- 
ticles and photons with the satellite, possibly 
complemented by exposures of the detector 
and associated structures to real particles and 
high energy photons, can be used to study so- 
lutions that will minimize the background. 

If the experiment is properly designed then 
the instrumental background will be low and 
the cosmic background important. Above ^ 
1 keV the dominant contributor to the cos- 
mic X-ray background is the extragalactic 
background associated to unresolved sources, 
mostly AGN. Sufficiently high spatial resolu- 
tion allows to resolve out a sizeable fraction of 
the sources producing the X-ray background 
(see Fig. |4j5). With a resolution of 5 arcsec 
(Half-Power-Ratio, HPR) it is possible to re- 
solve out about 80% of the background, pro- 



vided of course sufficient counts are available 
to detect the sources. It should be noted that be- 
yond an angular resolution of 15 arcsec the re- 
solved fraction is not very sensitive to the res- 
olution, see Fig.|4]3. Another important point is 
that, to fully exploit the advantage of a large 
field of view, it is necessary that the high spa- 
tial resolution be available over the full FOV, 
polynomial optics (Burrows et al. 1992) can 
provide this important feature. Another impor- 
tant contributor to the background is the so 
called straylight, this is associated to X-ray 
photons from outside the field of view which 
end up in the focal plane after reflecting only 
once on the mirrors. The effect of straylight can 
be significantly mitigated by introducing a pre- 
collimator in front of the telescope as was done 
in the case of the XMM-Newton optics. 
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Fig. 5. Simulation of source and background spectra for a typical region at R200. Tne instrumental back- 
ground estimates come from current missions. The cosmic X-ray background (CXB) is modeled with 3 
components, for the soft X-ray background (SXRB) we adopt the modeling of the SXRB from McCammon 
et al. (2002), who have carried out the highest spatial resolution observation of the SXRB with a sounding 
rocket flight. The SXRB is modeled by 2 thermal components with temperatures of 0. 1 keV and 0.225 
keV both with solar abundances, normalizations come from Table 3 of McCammon et al. (2002). For the 
extragalactic background, comprising mostly unresolved AGNs, we assume a power-law of slope 1.4 and 
intensity 1/4 of that derived by De Luca Sl Molendi (2004), thereby assuming that 3/4 of the sources will be 
resolved out. The source surface brightness is assumed to be 3 x 10"" erg cm^^s^'arcmin"^, a typical value 
for cluster outskirts (see Sect. l2.1l for a detailed discussion), the extraction region is 100 arcmin^, a value of 
2 X 10-''cm"- is assumed for the equivalent hydrogen column density, the cluster emission is modeled by a 
thermal plasma with kT=3 keV and metal abundance Z=0.15 Z©, the exposure time was set to 100 ks. The 
model spectrum was convolved with response files (effective are and redistribution matrix) provided by the 
VKFXrteam. 



3.2. Observing and data analysis 
strategies 

An experiment design like the one described 
above contributes significantly in improving 
the sensitivity to low surface brightness emis- 
sion, however further steps need to be taken 
to reach cluster outer regions. This is quite 
apparent when looking at the spectral simula- 
tion reported in Fig. |5] (for details see the fig- 
ure caption). As can be seen background com- 
ponents of one kind or another dominate the 
spectrum at all energies. In the 1-3 keV range 
the source intensity is about 1/3 of the total, 
below 1 keV the galactic foreground domi- 
nates, while above 3 keV the residual extra- 



galactic and the instrumental background do. 
These are of course estimates, for real clus- 
ters things may be a little different, however 
we will inevitably have a background that out- 
shines the source. These are atypical condi- 
tions with respect to previous X-ray imaging 
missions. To make reliable measures will re- 
quire devising specific observing and analysis 
strategies. Clearly the strongest requirement is 
that the background be characterized as well 
as possible, ideally one would like to mea- 
sure the background associated to the source 
without the source, which is of course impos- 
sible. Considering that the instrumental com- 
ponent varies with time and that the galactic 
foreground varies with position on the sky, it 
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is important to observe the background almost 
at the same time and almost at the same loca- 
tion of the source. A similar strategy has been 
adopted, albeit for reasons different from the 
ones considered here, by the SWIFT XRT ex- 
periment. During each 1.5 hour orbit, SWIFT 
observes a source field and 3 or 4 background 
fields. Thus background fields are observed al- 
most simultaneously with the source field and 
with the same instrument set up. Moretti et al. 
(2010) have shown that under these conditions 
the instrumental background can be character- 
ized to the 3% level. Conversely, when back- 
ground fields from dififerent epochs are used, 
only a 10-15% level is achieved. The optimal 
solution that may be applied in a future mis- 
sion, or on SWIFT for that matter, would be 
to use as part of the background fields, sky re- 
gions close to the source and dark earth ob- 
servations. The former would allow to perform 
a spatial characterization of the galactic fore- 
ground, while the latter would permit a clean 
measurement of the instrumental background. 
Observations of both source and background 
fields need to be conducted to a high preci- 
sion. Relative systematic errors on the spec- 
tra need to be kept at the few percent level. 
This is not a trivial requirement to meet, par- 
ticularly since at this level of precision each 
detector element has to be considered as an 
independent detector. Assuming that each de- 
tector element will be calibrated to a relative 
precision of ^ 5%, systematics can be reduced 
to the desired level by viewing each sky ele- 
ment with a large number of detector elements. 
Observing strategies such as this have been 
used for decades in other bands of the elec- 
tromagnetic spectrum when the source signal 
is smaller than the background. As examples, 
one may consider ground based infrared obser- 
vations or cosmic microwave background mea- 
surements. 

The comparison of the source plus back- 
ground spectrum with the background spec- 
trum is typically done via subtraction. In recent 
years, workers concentrating on cluster outer 
regions (e.g. Snowden et al. 2008, Leccardi 
& Molendi 2008) are finding that modeling 
is more efifective. This is readily understood 
if one considers that the background is made 



of different components each capable of vary- 
ing independently of the others. Unless there 
are good reasons to believe that the particu- 
lar combination of background components as- 
sociated to the source and background fields 
are next to identical, it is preferable to model 
the different components allowing for varia- 
tions in relative intensity. Another issue that 
should be considered is that, under the atyp- 
ical conditions of cluster outer regions, the 
standard maximum likelihood estimators com- 
monly employed to derive physical parameters 
such as emission measure and temperature do 
not always work properly. In a recent paper 
(Leccardi & Molendi 2007), we have shown 
that the presence of a significant background 
component can lead to a substantially biased 
measure of the temperature. In the same paper, 
we describe a few quick fixes. Unfortunately, 
a general solution, based on a more powerful 
statistical estimator, has yet to be found. 

3.3. A budget for systematics 

Assuming that the above guidelines are fol- 
lowed, we expect to be able to maintain sys- 
tematic errors to within a few percent. In the 
following, we provide a breakdown of the ex- 
pected errors. A constant monitoring of the 
instrumental background by using the part of 
the detector not exposed to the sky plus dark 
earth and background field observations en- 
twined with source observations should al- 
low us to constrain this component to about 
1 % (as extrapolated from the results obtained 
on SWIFT XRT in Moretti et al. 2010). The 
extragalactic component of the cosmic back- 
ground is a residual component, comprising 
unresolved sources and possibly a diffuse com- 
ponent. For a typical flux limit of lO""" erg 
cm^^s"' in the 0.5-2.0 keV band, montecarlo 
simulations show that the cosmic variance for 
a 100 arcmin^ field is less than 1% of the resid- 
ual background component. The galactic fore- 
ground will be monitored by performing ob- 
servations of fields contiguous to the source 
field. Moreover, observations over several 100 
arcmin^ should allow us to characterize this 
component to about 3-5%. Finally, assuming a 
typical relative calibration accuracy of 5% on 
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individual detector elements and the applica- 
tion of substantial dithering, we expect to reach 
an overall relative spectral calibration of about 
1%. 



3.4. Detailed predictions 

Our goal is to resolve the physical properties 
of the ICM in the virial regions making proper 
use of the WFXT (FOV with Rwf ~ 30'). 

Our strategy is to define a set of obser- 
vations with reasonable exposure time (< 50 
ksec) that can allow the study of the virial re- 
gions through the spatial and spectral analysis 
with WFXT. 

First, we select objects with known X-ray 
properties (flux, temperature, dynamical sta- 
tus) that can be good candidates for a single 
WFXT exposure, i.e. with an expected R200 < 
Rwp - 30'. We can also relax a bit this assump- 
tion requiring however that a given exposure 
minimizes the risks in term of (i) problems of 
intercalibration with other X-ray observatories 
for measurements in known X-ray emitting re- 
gions, (ii) weak constraints on the X-ray prop- 
erties at /?20o due to the effect of unexpected 
large scale structures. 

We estimate 7?2oo from a given spectro- 
scopic measurement of the gas temperature by 
using the best-fit results in Arnaud et al. (2005, 
cf. Table 2; similar results in Vikhlinin et al. 
2006): 



7?2oo = 1714x(rgas/5keV)"-^£^ 



:'/27(!kpc (1) 



r 1 nO-5 

with E. - [Om(l + z) + i^Aj and Q.^ = I - 

Qa = 0.3. 

By applying the criterion 7?2oo < ^wf = 
30', we select 23 out of the 45 objects present 
in the flux-Umited sample of the brightest clus- 
ters in Mohr et al. (1999; see Fig.|6]l). More ob- 
jects can be included if off-axis exposures are 
considered, as requested for the Perseus cluster 
with a R200 of ~ 88 arcmin. 

The response matrix used for our simula- 
tions is obtained by convolving the redistribu- 
tion matrix with a mean eff'ective area, A{E), 
constructed by averaging the vignetting over 



o 

Q. 



cr 




0.05 



0.10 0.15 
redshift 



0.20 



Fig. 6. Predicted R200 for the 45 objects in Mohr et 
al. (1999) compared to Rwf = 30' (red solid line) 
and the expected ^200 for a typical cluster with a 
temperature of 4 and 8 keV (dashed lines). R200 are 
estimated from eq.[T] 



the whole field of view, i.e. 



AiE) = A„(E) 



2KJ^"""6d0V(E,0) 



max 



where Ao(E) is the energy dependent on-axis 
effective area, 6,„ax is maximum off-axis an- 
gle and V{E, 0) is the energy and off-axis de- 
pendent vignetting. The redistribution matrix, 
the on-axis effective area and the vignetting 
were kindly provided by the WFXT team. We 
use our own script with the response matrix to 
simulate a (source+background) and a (back- 
ground only) spectrum including in the latter 
one (i) a 1 per cent random fluctuation in ab- 
sorbing «H value and in the normalization of 
the instrumental background; (ii) a 5 per cent 
random fluctuation propagated to the normal- 
ization and temperature values of the two lo- 
cal background component (one absorbed, the 
other not), to the normalization and photon- 
index value of the CXB. The photon-index is 
allowed to vary between 1.4 and 1.6. We as- 
sume that 80 per cent of the CXB is resolved. 
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Table 1. Relative errors (in percentage) and deviations e from the input values at 90% confidence level 
on the parameters of interest (normalization K, plasma temperature T and metal abundance Z of an apec 
component in XSPEC -Arnaud 1996) after joint-fit analysis of spectra simulated with an exposure time of 
50 ksec. All the relative errors can be rescaled to different exposure times as ~ y/exp- Pio indicates ayS value 
increased by 20 per cent. CC (nCC) indicates a (no) Cooling-Cores Cluster. The fluxes / are in units of 
10"'' erg/s/cm~ in the band (0.1 - 2.4) keV and are collected from|http://bax.ast.obs-mip.fr/| 



ij 



inputs 



K 



fixed Z 



K 



Perseus (TURBOLENT/CC; 

r = 3.16,Z = 0.15 

T = 2,Z = 0.15 

j82o,r = 3.16,Z = 0.15 

;e2o,7' = 2,Z = 0.15 

A1689 (MERGING/nCC; 

7' = 5.05,Z = 0.15 

r = 2,Z = 0.15 

j82o,7^ = 5.05,Z = 0.15 

/320,7' = 2,Z = 0.15 



z = 0.0178,/ = 1137.3, r = 6.3keV,nH = 1.5s21;i?200 = 1.9Mpc = 88.2') 
8(+0.2e) 15(+0.5e) 16 (+0.36) 17 (-1.36) >100 

8 (-1.06) 8 (+0.76) 17 (-0.16) 13 (-0.46) 44 (-0.26) 

27(+l.l6) 38(-0.7e) 66 (+0.46) 55 (-1.46) >100 

25 (+1.36) 30(+0.3e) 51 (+0.96) 47 (-0.26) 5(-1.06) 

z = 0.1810,/= 14.5, r= 10.1keV,«H = 1.8s20;i?200 = 2.2Mpc = 12.2') 
6 (-0.76) 23 (+0.56) 20 (-0.96) 27(+l.l6) >100 

6(+0.l6) 7(<0.l6) 13 (+0.26) 6 (-3.26) 35 (-1.36) 

>100 >100 >100 >100 >100 

55 (+0.26) 14 (-3.76) >100 >100 >100 



The spectra are integrated for 50 ksec over 
an area of 100 arcmin^ and then jointly fitted 
in the range 0.3-6 keV. 

The surface brightness in the band (0.5 - 2) 
keV are obtained from the best-fit values in 
Table 2 of Mohr et al. (1999) by evaluating 
the model prediction at 7?2oo as estimated in 
equation 1. A more conservative estimate of 
the surface brightness is obtained by increasing 
the p value by 20 per cent, faking an expected 
steepening of the surface brightness profile in 
the cluster outskirts, as recent observations and 
simulations suggest (see Section 2). This cor- 
rection reduces the predicted surface bright- 
ness by a factor of 7 on average. 

All the simulated spectra assume a metal- 
licity of 0.15Zq and a temperature equal to 0.5 
(see Roncarelli et al. 2006) times the quoted 
value in Table 1 of Mohr et al. (1999). We 
also consider the cases with metallicity equal 
to O.O5Z0 and temperature of about 0.25 times 
the quoted values (i.e. between 1 and 2 keV). 

Our simulated spectra (e.g. Fig. [S]) show 
that we can reach typical uncertainties (90% 
level of confidence) of < 20% on the normal- 
ization K and temperature T of the thermal 



spectra (see Tab.[T]). Reasonable constraints (~ 
40%) on the metallicity Z can be obtained in 
the case the surface brightness profile in the 
outskirts is still well reproduced from the mod- 
els fitted to ROSAT¥S?C data. 

A steepening of the surface brightness pro- 
files, as expected from the work discussed in 
Sect. 2 and modeled here by increasing the 
value of the outer slope/? by 20 %, reduces sig- 
nificantly the level of accuracy to which we can 
constrain the physical parameters: about 60 per 
cent (relative error at 90% level of confidence) 
on K, 40 per cent on T, no constraints on Z. 

4. Future missions & WFXT 

In this section we provide an overview of mis- 
sions under study or construction that may pro- 
vide important contributions to the characteri- 
zation of cluster outer regions. There are 3 such 
missions namely SRG, XENIA and WFXT. The 
eROSITA experiment (Predehl et al. 2007) on 
board the Russian Spektrum Roentgen Gamma 
(SRG) satellite comprises 7 telescopes with a 
total on-axis effective area of 2000 cm^, an on- 
axis angular resolution of 25 arcsec and will 
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operate from an L2 orbit. XENIA (Hartmann 
et al. 2009) carries an X-ray imager and spec- 
trometer that would both be useful in char- 
acterizing cluster outskirts: the imager has an 
on-axis effective area of 600 cm^ and an on- 
axis angular resolution of 15 arcsec; the spec- 
trometer has an unprecedented spectral resolu- 
tion of a few eV, an on-axis effective area of 
about 1000 cm^ and an angular resolution that 
is limited by the pixel size of a few arcmin. 
WFXr (Murray et al. 2010) which, like XENIA, 
has been submitted to the Astro2010: The 
Astronomy and Astrophysics Decadal Survey, 
carries an X-ray imager comprising 3 tele- 
scopes for a total on-axis effective area of 6000 
cm^, and an on-axis angular resolution of 5 
arcsec (requirement for the half-power-radius 
of the PSF). A low earth equatorial orbit is 
forseen for both XENIA and WFXT. 

Both the XENIA and WFXT imager have 
two considerable advantages over eROSlTA, 
namely the low earth over the L2 orbit and the 
polynomial optics, which will result in a sub- 
stantial reduction of the instrumental and cos- 
mic X-ray background, respectively. In partic- 
ular, the WFXr imager will provide the charac- 
terization of the cluster outer regions in about 
1/10 of the time requested from XENIA, and 
will benefit from higher angular resolution. 
XENIA however, is in the unique position to 
complement the imager data with high spec- 
tral resolution data for relatively bright clus- 
ters. While eROSlTA is scheduled for launch 
in 2012, XENIA and WFXTars both at an early 
stage of development and have to be consid- 
ered as the next generation satellites for clus- 
ters studies. 



5. Summary 

Past and current X-ray mission allow us to ob- 
serve only a fraction of the volume occupied by 
the ICM. Indeed, typical measures of the sur- 
face brightness, temperature and metal abun- 
dance extend out to a fraction of the virial ra- 
dius. The coming into operation of the second 
generation of medium energy X-ray telescopes 
at the turn of the millennium, has resulted in 
relatively modest improvements in our ability 
to characterize cluster outskirts. Even though 



recent results from Suzaku show some im- 
provement, the most sensitive instrument to 
low surface brightness to have flown thus far 
is quite possibly the SWIFT XKI which, ironi- 
cally, never had cluster outer regions as one of 
its top scientific objectives. 

The construction of an experiment capa- 
ble of making measures out to /?200 is well 
within the reach of currently available technol- 
ogy. What is required is an experiment design 
that will minimize the background, both instru- 
mental and cosmic, and maximizes the grasp, 
i.e. the product of effective area and FOV. 
Since cluster emission in the outskirts will be 
background dominated, instrument design and 
observational strategy should also allow for a 
meticulous characterization of the background. 
Detailed simulations based on realistic esti- 
mates of the different spectral components and 
of the precision with which the may be de- 
termined shows that an experiments such as 
the one we envisage will allow a solid char- 
acterization of cluster outskirts. From what we 
can surmise WFXT is, already designed to meet 
most of the requirements which are necessary 
to characterize cluster outskirts, and should 
have no major difficulty in accommodating the 
remaining few. 
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